Optomechanical sideband cooling of a micromechanical oscillator close to the 

quantum ground state 



R. Riviere/'QS. Deleglise*, 1 - 2 S. Weis*, 1 - 2 E. Gavartin, 2 0. Arcizet, 3 A. Schliesser, 1 ' 2 and T.J. Kippenberg 1 ' 2 <[j] 

1 Max Planck Institut fur Quantenoptik, 85748 Garching, Germany 
2 Ecole Polytechnique Federate de Lausanne (EPFL), 1015 Lausanne, Switzerland 
s Institut Neel, 38042 Grenoble, France 

Cooling a mesoscopic mechanical oscillator to its quantum ground state is elementary for the 
preparation and control of low entropy quantum states of large scale objects. Here, we pre-cool a 
70-MHz micromechanical silica oscillator to an occupancy below 200 quanta by thermalizing it with 
a 600-mK cold 3 He gas. Two-level system induced damping via structural defect states is shown 
to be strongly reduced, and simultaneously serves as novel thermometry method to independently 
quantify excess heating due to a cooling laser. We demonstrate that dynamical backaction sideband 
cooling can reduce the average occupancy to 9 ± 1 quanta, implying that the mechanical oscillator 
can be found (10 ± 1)% of the time in its quantum ground state. 



PACS numbers: 42.65. Sf, 42.50.Vk 



The quantum regime of mechanical systems has re- 
ceived significant interest over the past decade [TH3] . Me- 
chanical systems cooled to the quantum ground state 
may allow probing quantum mechanical phenomena on 
an unprecedentedly large scale, could enable quantum 
state preparation of mechanical systems and have been 
proposed as an interface between photons and station- 
ary qubits. To achieve ground state cooling, two chal- 
lenges have to be met: first, most mechanical oscillators 
have vibrational frequencies f2 m /27r < 100 MHz, such 
that low mode temperatures T c g are required to achieve 
Ml m > fceTeff (H is the reduced Planck constant and fee 
the Boltzman constant). Second, quantum limited mea- 
surements of mechanical motion must be performed at 
the level of the zero point motion, x zp f — yJh/2m c sQ, m 
in order to probe the state of the oscillator of mass m e s- 

Recently, a piezomechanical oscillator has been cooled 
to the quantum regime 0. Due to its GHz resonance 
frequency, conventional cryogenics could be employed for 
cooling, while it was probed using its piezoelectrical cou- 
pling to a superconducting qubit. In contrast, cooling 
schemes based on radiation pressure dynamical backac- 
tion as proposed [5J [7] and recently demonstrated [8HTU] 
can be applied to a much wider class of nano- and mi- 
cromechanical oscillators. This optomechanical scheme 
is based on parametric coupling of an optical and me- 
chanical resonance and simultaneously allows sensitive 
detection of mechanical motion. In analogy to the case of 
trapped ions [TT] , dynamical backaction sideband cooling 
[l2Tfl"5] can be used to reach the quantum ground state. 

Despite major progress, ground state cooling using 
this approach has remained challenging, owing to insuf- 
ficiently low starting temperatures or excess heating in 
the optical domain [1618 , while microwave experiments 
have been impeded by the residual thermal occupancy in 
the microwave cooling tone 19 or weak optomechani- 
cal coupling [20]. Here we demonstrate an experimental 
optomechanical setting that solves these challenges. 




FIG. 1: Cooling a micromechanical oscillator, (a) High Q 
mechanical and optical modes are co-located in a silica mi- 
crotoroid. The simulated displacement pattern of the me- 
chanical radial breathing mode (RBM) is shown; the opti- 
cal whispering gallery mode (WGM) is confined to the rim. 
(b) Thermalization of the RBM to the temperature of the 
3 He gas, with the lowest achieved temperature corresponding 
to an occupancy of the RBM below 200 quanta, (c) Opti- 
cal setup used for displacement monitoring of the mechanical 
mode, based on homodyne analysis of the light re-emerging 
from the toroid's WGM (see text for detailed description). 



We use silica toroidal resonators, which support whis- 
pering gallery modes (WGM) of ultrahigh finesse co- 
located with a low loss mechanical radial breathing mode 
(RBM) [211 [22] and large mutual optomechanical cou- 
pling. The devices used here have been optimized for 
narrow optical linewidths k, and moderately high me- 
chanical frequencies O m , thereby operating deeply in the 
resolved sideband regime (f2 m sa 27r • 70 MHz > 10k), 
while at the same time the pillar geometry was engineered 
for low mechanical dissipation [T51[23] (Fig. 1). 

For the cryogenic laser cooling experiments, we sub- 
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ject these samples directly to a 3 He gas evaporated from 
a reservoir of liquid 3 He recondensed before each experi- 
mental run. At a pressure of ~ 0.7 mbar, the gas provides 
a thermal bath at a temperature of ca. 600 mK. However, 
it is essential to verify thermalization of the toroid to the 
exchange gas. To this end, a low-noise fiber laser (wave- 
length A ps 1550 nm) is coupled to a WGM using a fiber 
taper positioned in the near field of the mode via piezo- 
electric actuators (Attocube GmbH). Using techniques 
described previously PH 023 , the displacement fluctua- 
tions of the RBM can be extracted and used to infer its 
noise temperature. As shown in Fig. [I]}), it follows the 
temperature of the helium gas down to temperatures of 
600 mK for weak probing (i.e. < 1 /jW input laser power). 

For the optomechanical sideband cooling, we employed 
a frequency-stabilized Ti:sapphire laser (A ps 780 nm), 
and a homodyne detection scheme 16J for quantum- 
limited detection of mechanical displacement fluctuations 
(Fig.Jlp). Importantly, for the high Fourier frequencies of 
interest the Tksapphirc laser is quantum limited in am- 
plitude and phase. Spectral analysis of this signal pro- 
vides direct access to the mechanical displacement spec- 
trum, from which the mechanical damping and resonance 
frequency can be derived. The spectra are calibrated in 
absolute terms j!6J [25] by applying a known frequency 
modulation at a fixed frequency close to the mechanical 
resonance frequency to the laser using an electro-optic 
modulator (EOM). After the acquisition of each spec- 
trum, the detuning of the laser from the cavity resonance 
is determined by sweeping the modulation frequency of 
the EOM and recording the demodulated homodyne sig- 
nal with the network analyzer |26j . 

Before studying radiation-pressure induced effects, we 
have carefully analyzed the influence of the bath temper- 
ature on the RBM's properties. The vitreous nature of 
silica leads to a strong temperature dependence due to 
the presence of structural defects modeled as two-level- 
systems (TLS) [H [27]. Relaxation of the TLS under 
excitation from an acoustic wave modifies the complex 
mechanical susceptibility, leading to a change in mechan- 
ical resonance frequency J7 m and a change of the damping 
rate T m = £l m /Q m with mechanical quality factor Q m . 

Two different relaxation regimes have to be considered 
[2"B] for sample temperatures T between 0.6 K to 3 K: 
tunneling-assisted relaxation [29J, [30] , and single phonon 
resonant interaction 30J . Thermally activated relaxation 
[31) dominates the frequency shift at temperatures above 
3 K, but is negligible in the temperature range at which 
the laser cooling experiments arc performed. In the pres- 
ence of tunneling relaxation ("tun") and resonant inter- 
actions ("res") the mechanical oscillator properties can 
be expressed as 

On(T) = fl m + SQ tun (T) + 5il Tes (T) (1) 
r m (T)/a„ = Q m \T) = Q2 + QtunCH + Q7ol(T), (2) 
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FIG. 2: TLS-induced change of the resonance frequency (fi m ) 
(a) and inverse mechanical quality factor (Q^, 1 ) (b) of the ra- 
dial breathing mode. Measured data (red points) agree well 
with the models (red lines) described in eqs. Q and (2). 
Subtraction of the temperature-independent clamping damp- 
ing (line iii) yields the theoretically possible material-limited 
damping values (blue points and line ii). At very low temper- 
atures, damping by resonant interaction with TLS (iv) would 
be dominant. The model parameters are given in ref. |28| . 



where ^mQ^a * s the damping rate due to the clamp- 
ing of the resonator to the substrate, dominating r m at 
room temperature. The respective temperature depen- 
dencies in the relevant regimes of TLS damping are de- 
tailed in [25]. For the lowest temperatures of 600 mK, Q nl 
reaches ~ 10 4 sufficient to enable ground state cooling 
since Q m /ni ^> 1 and fhiQ m /Q m <C K (nj is the initial oc- 
cupancy) |32j . Moreover, the well- understood tempera- 
ture dependence of the TLS-induced effects [Eqs. ([T])-([2|] 
enables its use as a "thermometer" of the sample temper- 
ature T after a calibration measurement as shown in Fig. 
[2] has been performed once. Importantly, this method 
can reveal excess heating independent of the RBM's noise 
temperature. 

We next studied optomechanical cooling [BTfTU] by per- 
forming a series of experiments in which mechanical dis- 
placement noise spectra were recorded while varying the 
laser detuning A = u>\ — uj c . Here, w\ is the laser's (angu- 
lar) frequency and ui c the WGM frequency, taking tem- 
perature and static radiation-pressure induced shifts into 
account. The mechanical mode's frequency and damp- 
ing shows a strong detuning dependence (Fig. [3]) owing 
to the dynamic in-phase and quadrature response of the 
radiation pressure force with respect to the mechanical 
motion. 

To accurately model radiation-pressure induced dy- 
namical backaction [6] for the present microresonators, 
we have to additionally take into account that backscat- 
tering of light can couple WGMs with opposite cir- 
culation sense [551 1331 EH]- The mutual coupling of 
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the clockwise (a cw ) and counterclockwise (a ccw ) orbit- 
ing modes lifts the degeneracy leading to new eigen- 
modes of the system, i.e. a + = (a ccw + a cw )/V2 and 
a_ = (a ccw — a cw )/-\/2, where the new eigenfrequencies 
lu± = ui c =F 7/2 are split by the mutual coupling rate 7. 
During a detuning series as reported here, both modes 
are populated by the driving field s- m with a mean field 
o± = \J k cx /2 L±(A)s- m , where Pj n = |si n | 2 fex;i is the 
driving laser power, k cx the coupling rate to the fiber ta- 
per, |a-t| 2 the mean photon population in the new eigen- 

modes and L±(A) = (-»(A ± 7/2) + k/2) _1 the modes' 
Lorentzian response. 

In the context of cavity optomechanics, it is important 
to realize that three-mode interactions [35] can be ne- 
glected in the present configuration j5S]. The radiation 
pressure forces induced by the light in these modes can 
therefore simply be added. The usual linearization pro- 
cedure [3(3 [37] then yields an inverse effective mechanical 
susceptibility of 

(xeff my 1 = TOeff -n 2 - ir m n - io m /(0)) (3) 

modified by dynamical backaction according to 
/(ft) = 2g 2 £ K\ 2 {L a {K + n) - (i CT (A - fi))*) (4) 

cr=± 

with the vacuum optomechanical coupling rate |25j go = 
Gx zp { and G = dcj c /da;. For moderate driving powers 
[32] , the susceptibility of the mechanical mode is the one 
of a harmonic oscillator with effective damping rate and 
resonance frequency of 

r eS » r m (T) + Re [/(n m )] (5) 
n cS ^n B1 (T) + im[f(n m )}/2. (6) 

For the samples studied in the following, a coupling 
rate of |<?n| ~ 27r x (1.2 ± 0.2) kHz is determined from the 
coupling parameter \G\ = uj c /R s» 2ir x 16GHz/nm and 
effective masses m e g = 20 ± 5ng. Figure [3] shows the re- 
sults of a detuning series, which was taken with an input 
laser power of 2 mW, with the temperature of the 3 He gas 
stabilized to T cryo = 850 mK at a pressure of 2.8mbar. 
The excellent stability of both the laser and the cryogenic 
coupling setup allowed us to perform the series without 
active stabilization of the laser detuning A and the cou- 
pling k cx , which is determined by the sub-micrometer 
gap between the coupling taper and the toroid. 

The coupled fit of the data using the model of eqs. 
([3|-([6]) is enabled by the precise calibration of the laser 
detuning by sweeping the modulation frequency . We 
choose to adjust the parameters of the fit primarily (rela- 
tive weight 0.9) to the optical spring effect, since the me- 
chanical resonance frequency can be extracted from the 
spectra with higher accuracy than the damping rate. The 
obtained fit parameters k, 7, and s; n are found to be in 
good agreement with independent results deduced from 
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FIG. 3: Effective resonance frequency (a) and linewidth (b) 
of the RBM when a 2 mW-power laser is tuned through the 
lower mechanical sideband of the split optical mode (inset). 
Blue points are measured data extracted from the recorded 
spectra of thermally induced mechanical displacement fluctu- 
ations, solid lines are a coupled fit using the model of eqs. 
(|3|-l|6|, taking into account heating of the cavity due to ab- 
sorbed stray and intracavity light, modifying the mechanical 
properties via the temperature dependence of the TLS. 

the frequency modulation measurement (kw 27t x 6 MHz, 
7 2n x 30 MHz) and the measured laser power. 

The excellent quality of the fit, together with the mea- 
sured temperature dependence of the TLS effects on the 
mechanical mode, furthermore allows us to extract the 
temperature T of the sample. Importantly, for large de- 
tunings |A| 3> k, the TLS thermometer reveals an in- 
crease of the sample temperature by AT stray « 220 mK 
corresponding to Q m (T) = 5970, which we attribute to 
heating induced by absorbed stray light scattered from 
defects on the fiber taper, which were observed to aggre- 
gate upon its production. 

As the laser is tuned closer to resonance, more light is 
coupled into the WGM and 

T fa T cryo + AT stray + AT W gm > (7) 

where ATwgm = /3K a bs(|a+| 2 + \a-\ 2 )huj\ denotes the 
increase in temperature following the cavity's doublc- 
Lorentizan absorption profile, K a b s ^ K — k ox is the 
photon absorption rate and j3 the temperature increase 
per absorbed power. Operating deeply in the re- 
solved sideband regime [15], only little optical power 
(~ K a b s |a + | 2 tvjjy < P in /1300) can be absorbed in the cav- 
ity when A = — f2 m — 7/2, leading to a modest tempera- 
ture increase of ATwgm ~ 70 mK. 
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Importantly, we can test the consistency of the de- 
rived detuning-dependent quantities — T, T m (T), and 
r e ff — by comparing the expected effective temperature 
of the RBM due to optomechanical cooling, i.e. T e g = 
T ■ r m (T)/r e fj, with the effective temperature derived 
from noise thermometry via integration of the calibrated 
noise spectra |16) . Figure 4a) shows this comparison for 
the detuning series discussed above. Using the model of 
cqs. ([2])- ^ adjusted to the data of Fig.[3j we obtain good 
agreement for the effective temperatures obtained in both 
ways. To achieve this level of agreement, it is necessary 
to take into account the optomechanical de-amplification 
of the laser phase modulation used for calibrating the 
mechanical fluctuation spectra in absolute terms: as was 
shown in a recent study [38] , the transduction of a phase 
modulation of depth Sip at a frequency f2 mo d/27r is mod- 
ified by a factor |x e ff(^mod)/Xm(^mod)| in the presence 
of dynamical backaction, where x m (f2) is the bare me- 
chanical susceptibility. 

Figure 4b) shows the same comparison for a cooling 
run at a high laser power (4 mW) , for which we ob- 
serve slightly increased heating by AT stra y ~ 400 mK, 
while additional heating by AXwgm could not be dis- 
cerned in this measurement. In spite of the reduced 
mechanical quality factor Q m (T cryo + AT s t ra y) = 4880, 
the lowest extracted occupancy is n = 10 according 
to the detuning series fit. The lowest inferred noise 
temperature of a single measurement is even slightly 
lower, corresponding to n = 9 ± I, where the uncer- 
tainty is dominated by systematic errors, which we es- 
timate from the deviations of the effective temperature 
derived in the two independent ways described above. 
Note that this occupancy implies already a probability 
of P(n = 0) = (1 + n) _1 = (10 ± 1)% to find the oscilla- 
tor in its quantum ground state. From this measurement, 
we can also extract a total force noise spectral density of 
Sff — (8±2 fN/yTIz) 2 driving the oscillator to this occu- 
pancy, corresponding to a imprecision-backaction prod- 
uct [H H5] of y/ S ^ S FF = ( 49 ± 8 W 2 if the entire 
present force noise is (conservatively) considered as mea- 
surement backaction. Indeed, wc estimate about 40% of 
the force noise to originate from the Langevin force at 
the temperature of the cryostat, and 60% from the ex- 
cess Langevin force caused by the laser-induced heating 
by AT s t r ay EH- -Due to the resolved-sideband operation, 
force noise due to quantum backaction (as yet only ob- 
served on cold atomic gases [39]) is expected to be nearly 
two orders of magnitude weaker and therefore negligible. 

It is realistic to significantly reduce the occupancy by 
higher cooling powers and improvements of go/r m , which 
can be achieved by engineering mechanical modes [23J for 
smaller mass and lower damping. For occupancies n < I , 
we anticipate that individually resolved anti-Stokcs and 
Stokes sidebands [HI [T5] of an independent readout laser 
will display a measurable asymmetry of n/(n+ 1) arising 
from the non-zero commutator of the ladder operators 
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FIG. 4: Cooling factor (T cryo + AT s tray)/T e ff and phonon oc- 
cupancy number of the RBM versus laser detuning A for (a) 
Pin = 2mW and (b) Pj n = 4mW. Points are phonon occu- 
pancies derived from the measured noise temperature, while 
solid lines correspond to the occupancies expected from the 
dependency of sample temperature and (intrinsic and effec- 
tive) damping extracted from the detuning series as shown in 
Fig. [3j A minimum phonon number of n = 9 ± 1 is obtained. 
The inset shows a mechanical displacement noise (DN) spec- 
trum at the optimum detuning (A = — f2 m — 7/2), illustrating 
the signal-to-noise ratio achieved despite the low occupancy. 



describing the mechanical harmonic oscillator in quan- 
tum mechanical terms. 
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Supplementary information - Optomechanical sideband cooling of a micromechanical 

oscillator close to the quantum ground state 



TWO-LEVEL SYSTEMS 



Tunneling systems in SiC>2 play an important role in cryogenic operation of silica mechanical oscillators. They lead 
to a temperature dependent frequency shift (via a change in speed of sound) of the considered mechanical mode and 
temperature dependent mean free path of phonons in SiC>2 affecting the mechanical quality factor. We will discuss 
these effects in the following and give the relevant formulae which have been used to fit the data in the main part of 
the manuscript. An extensive study of TLS effects can be found in [2 [2]. 

As hrst considered by L. Pauling in 1930 [3 , tunneling of atoms occurs in solids with a certain degree of disorder, 
where in the local environment of an atomic site, several potential minima exist. This can be the case in the vicinity 
of defects in crystals or, more frequently, in amorphous materials. For amorphous solids at low temperatures, the 
tunneling dynamics can be well captured in a simple model consisting of an ensemble of two level systems (TLS) 



each of which is described by a generic double- well potential (Fig. SI). This potential is parametrized only by the 
barrier height V , the initial energy asymmetry Ai and the spatial separation between the two potential minima d. A 
tunneling coupling strength 

A = M7 e~ A (SI) 

with the intrinsic oscillation frequency fio within the individual atomic sites can then be deduced, with the tunneling 
parameter 



x~ 2mVd «so\ 
A "V^2 (S2) 

depending on the atomic mass m. Due to this tunnel coupling, the new eigenmodes of the coupled system exhibit an 
energy splitting of 



E=JAj + Al. (S3) 




FIG. SI: Double well potential with relevant levels and naming convention 

Phonons couple to TLS via their strain field that leads to a deformation of the TLS potential (notably leading to 
a change in the barrier height V). As a consequence, the TLS are driven out of thermal equilibrium and relaxation 
processes will exchange energy with the heat bath. Transitions between the two energy levels are induced by several 
distinct processes that become dominant in different temperature regimes. 
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• At very low temperatures (T < E/k&), the density of thermal phonons is low, such that relaxation processes 
do only play a minor role. Here a significant population imbalance between lower and excited state exists, and 
the most efficient transition mechanism is resonant absorption of phonons of a frequency fi ra = E/H. This 
mechanism shows — as in the case of other two level systems — a saturation behavior [3]. 

• At temperatures T > E/k-Q (typically a few Kelvin), the number of thermal phonons has increased to a level 
at which Raman processes involving a tunneling through the barrier become predominant. It is mainly this 
temperature range that will be relevant for the description of the phenomena seen in the present cooling exper- 
iment. 

• At even higher temperatures thermally activated relaxation dominates. In this case multi-phonon processes 
with an excitation across the barrier take place. 



Relaxation contribution 



If relaxation is the dominant process, the general expression for the mean free path of a phonon of frequency f2 m is 
given by [H |2] 

rl<T) =M(-S) 4 ^T^ (a '< A)da ' dA - (S4) 

The integration is performed on all TLS that can interact with the phonon. Here, P(A 1; A) is the volume density of 
TLS with energy asymmetry between Ai and Aj + dA x and tunnel parameter between A and A + dA, 

n ° = e E/k B T + i ( S5 ) 

is the thermal equilibrium Boltzmann repartition function, c s the speed of sound, p the mass density of the solid, 
r the relaxation time of the individual TLS and B the coefficient linking a deformation Se to a change of E via 
SE = 2B(Ax/E)Se. 

A mechanical quality factor of 

Q m 1 (r) = ^^ (S6) 

can then be deduced. For the corresponding relative change in the speed of sound (i.e. frequency shift of a mechanical 
resonance) one obtains from the Kramers-Kronig relations 

«WT) = -fe // {^^—^— P(A l5 A)dA ld A. (S7) 



2pc s 2 77 V 9E J £ 2 1 + ft^r 2 



Tunneling- assisted relaxation 



Within the framework of the so-called tunneling model p] [2] the relaxation time is given by 

E 2 



l A 



(S8) 



with the maximum relaxation rate 



^ <S9) 

Parametrizing the integrals in terms of the energy splitting E and the parameter u — r^ 1 /r^ 1 yields [HE] 

" (-H) [ OTI d " dE (S10) 



and 

m ^ {T)= -^rL \-^)Lu^^j u ^ (S11) 

where it is assumed that the density P(E, A) = P is constant, which is consistent with experiments. A prominent 
feature in this regime is a plateau of the quality factors for temperatures of a few Kelvins with Q values of 

Q~} . =1*-^. (S12) 
^plateau 2 pc 2 y ' 



Thermally activated relaxation 
At higher temperature the rate is given by the Arrhenius law and only depends on the energy barrier height, 

T -l - T -l p -V/k B T 
T th — T e I 

where tq represents the period of oscillation in individual wells p] H] . 



(S13) 



Resonant processes 

For resonant interaction between phonons and TLS, it can be shown that [HE] 



fcs(T) = ^tanhf^%) (S-Li, 



pc 2 \2k B T 
pel U {t 



( yn res ( T ) = ^_i n (_) ! (S15) 



where To is a reference temperature. While resonant processes do not significantly contribute to the mechanical 
quality factors in our experiment, the frequency shift is dominated by resonant processes. 

Fitting Parameters for Figure 2 

The curves shown in figure 2 of the main manuscript have been fitted with the equations given in the previous 
sections. For the frequency shift the sum of the tunneling relaxation and the resonant contribution has been taken into 
account. The latter dominates this effect up to about T — 2K. The contribution of thermally activated relaxation has 
been omitted since it doesn't contribute significantly in the considered temperature range. Fitting of the Q-dependency 
has been done using the sum of tunneling relaxation, resonant contribution and a constant offset accounting for the 
clamping losses (Q^a)> ^- e - ^ oss °^ acous tic energy due to leaking into the substrate for this particular toroid. Here 
the resonant contribution plays a minor role. 

For the curves shown in Fig. 2 of the main manuscript, we used the material parameters 

c s = 5800 m/s 
p= 2330 kg/m 3 , 



the measured resonance frequency 



as well as the adjusted parameters 



tt m = 2tt x 76.3 MHz, 



B = 1.1 x 10~ 19 J 
P Qm = 2.5 x 10 45 m~ 3 
P 0m = 4.6 x 10 45 m" 3 . 
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For the fitting of the two curves (mechanical quality factor, resonance frequency shift) two different values for P 
had to be used. Given that the two traces are governed by two different regimes, small differences in the density of 
contributing TLS to the two effects seem to be justified. The literature [6] value of the dimensionless parameter C — 
PB 2 /(pel) = 3.010~ 4 shows a reasonable agreement with the parameters of the resonance frequency (C^ m = 7.110 -4 ) 
and damping (Cq^ = 3.910~ 4 ) fits. 

DYNAMICAL BACKACTION IN THE PRESENCE OF MODE SPLITTING 

In the framework of coupled-mode theory [7] , the two coupled counterpropagating modes [H [5] in a WGM resonator 
can be described by the equations of motion (in a frame rotating at the laser frequency) 



a CC w(£) = («(A - Gx(t)) - K/2)a ccw {t) + i-a cw (t) + y/ryis- m (t) 
a cw {t) = (i(A - Gx(tj) - K/2)a cw (t) + i-a ccw (t). 



(S16) 
(S17) 



Here rj c describes the coupling parameters defined via rj c = - h + Kg , where k cx describes the output coupling rate, 
whereas kq denotes the intrinsic loss rate of the cavity. 
The fields in the system's new eigenmodes 



a+ = (flccw + a cw )/\/2 
a~ = (a C cw - a cw )/V2 



(S18) 
(S19) 



exert a radiation pressure force of 



F rp = -/iG(|a + (<)| 2 + M*)| 2 ), 



(S20) 



since the spatial shape of cross-term 2Re(a+(t)a*_(t)) has an azimuthal dependence oc cos(mip) sm(mip) (m is the 
angular mode number), averaging to zero when projected on the azimuthally symmetric RBM. The coupled optome- 
chanical equations of motion can therefore be written as 



,Mn = [ i [A-Gx(t) + ^)^) fl+ (() 



d_(i) = (i (A -<?*(*) -J)-!)a_(i) + 
™ cff (x(t) + T m x{t) + n 2 m x(t)) = -hG{\a+{t)\ 2 + \a-(t)\ 2 )+6F(t), 



2 

T] C K 




Sin(*) 
Sin(t) 



(S21) 

(S22) 
(S23) 



where SF(t) is an external force, e. g. the thermal Langevin force. 
We then apply the usual linearization 

a±(t) = a± + Sa±(t) 
x(t) = x + Sx(t) 



(S24) 
(S25) 



assuming \a±\ ^> \Sa±(t)\ and |x| ^> For the large mean occupancy of the modes and the mean mechanical 

displacement, we then obtain 



a + 



y/T) c K/2i 



-i(A + 7 /2) + K /2 

y/rj c K/2s iD 
-t(A - 7/2) + k/2 



^ VcK /2L + {A) 



HG 
m cB tt 2 , 



(S26) 

(S27) 
(S28) 
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The average displacement x induces a small static frequency shift, as does the (usually dominant) static shift due to 
absorption- induced heating [10) . which are both absorbed into the mean detuning A = uj\ — (u) c (T) + Gx). One then 
obtains the equations of motion of small fluctuations, 

6a+(t) = (»(A + !)- |) Sa+{t) - iGa+Sx(t) (S29) 
Sa-(t) = (i(&-T\- J) 6a-(t) - iGa^6x(t) (S30) 



27 2 

m Gff (Sx(t) + T m 6x{t) + O^fa(i)) = ~hG (a* + Sa+(t) + a+ (Sa+(t))* + a*_5a_(t) + a_ (<5a_(f))*) + 5F(i). (S31) 
Fourier transformation gives 

— z (A + 7/2 + fi) + K/2 

<5a_ (fi) = -iGa 5x{Sl) = _ ^ 

V 7 -i(A- 7 /2 + fi) +k/2 V y V ' y ' 

8x(n)/xm(fy = -KG (a* + Sa + (+Q) + a+ (Sa+(-Q))* + a*_Sa-(+Cl) + a_ (<5a_(-fi))*) + 5F(fi). (S34) 

With 

= 7 02 1 -nr r~cr2\ ( S35 ) 

m eff {-ll 2 - iilT m + fi£J 

Solving equations (S32 - S34) for 5x yields 

5x(Sl) = t > (S36) 

1/Xm(fi) - *fiG 2 (ja+| 2 LL+(A + fi) - (£+(A - SI))*) + |a_| 2 (i_(A + fi) - (£_(A - fi))* 1 

so that we can write 

— = — L- - iftG 2 (|a + | 2 (i + (A + fi) - (L+(A - fi))*) + |a_| 2 (i_(A + 17) - (£_(£ - fi))*)) (S37) 

and, in the regime of weak optomechanical coupling [llj 

T cff w T m + 2x 2 pf G 2 Re (|a+| 2 (i+(A + SI) - (£+(£ - fi))*) + |a_| 2 (i_(A + fi) - (£_(£ - fi))*)) (S38) 
fioff « fim + ^z P fG 2 Im (|a+| 2 (l+(A + fi) - (L+(A - fi))*) + |a_| 2 (i_(£ + fi) - (£_(£ - fi))*)) . (S39) 

CALCULATION OF THE IMPRECISION-BACKACTION PRODUCT 

In the context of quantum measurements [12j . it is interesting to characterize the sources of noise responsible for 
the mechanical displacement measurement uncertainty. For a given mechanical spectra, the measured (double-sided, 
symmetrized) spectral density of displacement fluctuations is given by 

5- oas (fi) = 5^(fi) + | Xcff (fi)| 2 SFF(fi) (S40) 

where S , "° p (fi) describes the measurement imprecision due to apparent displacement fluctuations, which are actually 
caused by noise in the displacement transducer itself. Sff(SI) is the force noise acting on the mechanical oscillator, 
and Xeff(^) its effective mechanical susceptibility. It is particularly interesting to evaluate these quantities for the 
lowest occupancy obtained at the optimum detuning of A = — fi m — J and at the Fourier frequency fi = fi m . 
In our experiment, the measurement imprecision is dominated by shot noise, and we extract a value of 

S™ p = 5™ p (fi m ) = (3.2 x lO-^m/x/Hz) 2 

from the fit to the background of the recorded mechanical spectrum (Fig. 4 from the main manuscript). Its measured 
linear dependence on the laser input power P ln shows that it is strongly dominated by the quantum noise of the 
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input laser. This behavior is indeed expected at the frequencies of interest in our work, where classical quadrature 
fluctuations are negligible in Tksapphire lasers. 

The thermal force noise (for j^- ^> 1) driving the mechanical oscillator is given by 

= S F h F c (f! m ) = 2m cS k B TT m (T) (S41) 

given by the fluctuation-dissipation theorem. In the presence of dynamical backaction, we can estimate this force 
noise from the more directly measured linewidth T c g and noise temperature T e ff using T c ff = T ■ r m (T)/T e ff, and 

S F h F c = 2m cff fc B T cff r cff (S42) 

It evaluates to 

S& e = (8 ± 2 x 10~ 15 N/Vlz) 2 

where T e ff and T e g are extracted from the fits to the detuning series, evaluated at the detuning A = — fi ro — ^ as 
described in the main manuscript. This value gives a conservative estimate of the classical measurement backaction, 
considering effectively all force noise present in the system (including thermal noise due to the non-zero cryostat 
temperature) as a classical backation of the measurement. 

A less conservative estimate on the backaction of the actual displacement measurement using the laser coupled to 
the WGM can be made by separating two different contributions in the force noise, 

cthe _ c cr y° i cba iQ<AO\ 
Opp — ^FF ' ^FF' yo^O) 

where S FF ° is the Langevin force noise due to the finite cryostat temperature T cryo and S FF the thermal backaction 
in the form of excess Langevin force noise due to the heating of the cavity by laser light. S FF gives an estimate 
of the classical perturbation of the system by the measurement, the classical excess backaction, which is technically 
avoidable. 

The thermal force noise originating from the bath 

SpjT = 2m cff A:BT cryo r m (T cryo ) (S44) 

is estimated to 



Sp r F yo = (5 ± 1 x 10- 15 N/VHz) 2 . 

T cryo and r m (T cryo ) are extracted from independent low input power measurements where the RBM is thermalized 
to the cryostat temperature. 

Consequently, the excess classical backaction evaluates to 

££1 = (6 ± 2 x 1CT 15 N/%/tfe) 2 

and accounts for 60 % of the thermal force fluctuations driving the mechanical oscillator. 

In addition to classical backaction, the quantum fluctuations of the intracavity photon number give rise to a quantum 
measurement backaction for which the force noise is given by 



D FF = ^FF I s 'm J 02 — ~=~Tj l 040 J 



in the case of high resolved sideband factor ^* 3> 1 [13] and at the detuning of interest. It is of the order of 
(1 x 10 _15 N/vHz) 2 in our case, negligible compared to the classical backaction. 

on \\~\T ( Tf\T - 

m„ 



Therefore, a conservative estimate of the imprecision-backaction product is given by (for ^> 1) 



S™ p (S F h F c + S£ F a ) w V^Sff = (49 ± 8)5 

In an ideal quantum measurement |12j . this product is equal to |, corresponding to the optimal compromise be- 
tween quantum imprecision and quantum backaction, both arising from the quantum fluctuations of the optical field 
quadratures. 
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As it is shown in the main manuscript, laser absorption heating, responsible for the classical excess backaction 
Spp, is mainly caused by scattered light off the tapered fiber being absorbed by the toroid (in our case by dust 
particles on the tapered fiber originating from particles in the air of our laboratory). It is thus within technical 
reach to strongly reduce this effect and perform measurements where light induced backaction would be dominated 
by quantum fluctuations alone. 
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